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P ■ Abstract 

Gravitinos with a mass in the keV range are an interesting candidate for warm 
dark matter. Recent measurements of the matter density of the universe and of 
^ ■ cosmic structures at the dwarf galaxy scale rule out the simplest gravitino models 
with thermal freeze-out. We construct a model where the decay of the messenger 
particles that transmit the supersymmetry breaking to the observable sector generates 
the required entropy to dilute the gravitino relic density by the required factor of a 
few to come in line with observations. The model is natural, and requires only that 
the coupling of the messenger sector to the standard model be set so that the decay 
happens at the appropriate time. 



1 Introduction 



We revisit the possibility that gravitinos could be the dominant component of non- 
relativistic matter in the universe. Especially interesting to us are scenarios where 
the gravitino is very light, with a mass in the keV range. Such a species would be 
"warm" dark matter. 

This possibility is interesting for two reasons, which we will now enumerate. The 
first is that such a light gravitino arises naturally in models of gauge-mediated su- 
persymmetry breaking M (GMSB) revived recently 0. They are stable, and thus 
could potentially pose problems for cosmology [§J |^, ||, which we will address in light 
of new more accurate measurements of both the matter density of the universe and 
on constraints on deviations from cold dark matter. 

Our second motivation is that recently cold dark matter has come under fire, and 
warm dark matter may in fact be favored. In particular, the predictions of the cold 
dark matter theory indicate that galaxies such as the Milky Way should have a large 
number of satellites, which are not observed. Warm dark matter could ameliorate 
this problem, as density fluctuations at the satellite scale and smaller are suppressed 
H |^|. In these models, the densities of halo cores are lowered, large halos have 
many fewer subhalos, and the subhalos that are present form in a top-down fashion, 
namely by fragmentation of larger clumps. This is contrary to the prediction of cold 
dark matter, where all relevant structures form in a bottom-up fashion by mergers 
of smaller structures ||. 

As data have improved, both in terms of the matter density of the universe, and 
in constraints on deviations from collisionless dark matter, the gravitino has silently 
stopped being a viable dark matter candidate as its relic density is too large in the 
simple thermal freeze-out scenario. We rescue the gravitino in GMSB models by 
invoking a natural model for entropy production, and thus provide again a candidate 
for warm dark matter, with a mass in the range 1.0 — 1.5 keV. 



2 Gravitino mass and relic density 

In GMSB models, the gravitino mass m 3 / 2 is related to the SUSY breaking scale 
Asusy as follows, 

m -^Rr)-°- 237 kev (w) 2 ' « 

where M p \ is the reduced Planck mass. We will be concerned with gravitino masses 
around 1 keV, thus SUSY breaking scales of 2 PeV, which are natural in GMSB 
models.[] 

■"■In this class of models, so-called messenger particles acquire a mass M and a supersymmetry- 
breaking scalar mass term F. The superparticles in the SUSY Standard Model acquire soft masses via 
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Gravitinos in this scenario interact fairly strongly, and thus freeze out quite late, 
at a temperature of order the weak scale or even colder. As they are of course highly 
relativistic, the computation of their relic density is relatively simple: all that needs 
to be done is to compute the effective degrees of freedom g* at freeze-out. Their relic 
density is then obtained, 

The value of g* at gravitino freeze-out has been determined to be in the range 90-140 
for a wide range of parameters 1|10|| . 

Finally we consider astrophysical constraints on the mass of the relic particle 
making up the dark matter. Particles in the eV range that undergo thermal freeze- 
out are ruled out as "hot" dark matter. The transition to "cold" dark matter occurs 
for thermal relics in the keV range. In fact, there is a lower bound on the mass of a 
warm relic DDI, 01, 



m 3/2 > 0.75 - 1.0 keV — ► Q 3/2 h 2 > 0.6, (3) 

based on the requirement that the small scale fluctuations are not in conflict with 
e.g. the Lya forest. This mass bound is valid for any thermal relic, but is modified 
for relics whose distribution functions deviate from equilibrium fl~2|l . 

We compare this bound on the relic density of gravitinos with current data on the 
density of the dark matter, which favors f^DM^ 2 ~ 0.15 [ 13f| , and in any circumstances 
Odm^ 2 < 0.3, indicating a serious discrepancy. One way to alleviate this is to allow 
for entropy production after gravitino freeze-out, which dilutes the gravitino relic 
density to acceptable levels. We discuss this possibility in the next section. 



3 Entropy production: simple model 

The production of entropy is required to dilute the relic density of gravitinos. As 
current data favor f^DM^ 2 ~ 0.15, and in any event Odm^ 2 < 0.3, we will consider a 
comoving entropy density increase of a factor of several. 

j|^2 J (-jg) , where g a are gauge coupling- 
constants and C a are quadratic Casimir. In the original model of gauge mediation, \f~F ~ M is 
suppressed relative to the primordial supersymmetry breaking scale Asusy by a loop-factor, and the 
gravitino mass is typically heavier than 50 keV But in models of direct gauge mediation ||, 
\fF ~ Asusy with no loop suppression. In particular, the model in Ref. ]7j achieves m 3 / 2 ~ 1 keV 
naturally. We thus can have standard weak-scale supersymmetry in the rest of the superparticle 
spectrum, with the caveat being that the gravitino is the LSP and all other superparticles decay to 
gravitinos and Standard Model particles. 
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In GMSB models there is a messenger sector around the SUSY breaking scale, 
which for keV gravitinos is of order 2 PeV. As our entropy production mechanism 
we will consider the out-of-equilibrium decays of the messengers via small couplings 
to the Standard Model. Assuming stability of the messengers, their relic density has 



been calculated fFl, and we simply quote the results here, 



n M h> ~ 10» f - Yu = ^ « 3.65 x HT" (JL) , (4) 

M VPeV/ s M VPeV/ ' V ; 

where M is the mass of the messenger particle, Ym is the frozen-out value of n/s, the 
ratio of number and entropy densities for the messenger particles, Sq is the entropy 
density today, and p c is the critical density today. 

In this section we will assume that the messengers decay instantaneously when 
the temperature of the thermal bath is To, and their Standard Model decay products 
thermalize instantaneously also. This approximation requires only the conservation 
of energy. Before decay, the energy density is 

p=^g*T i + 2 ^-g^MY M . (5) 

The first term is the radiation density and the second is the matter density of mes- 
sengers. This quantity is simply equated to the radiation density afterwords, with a 
higher temperature T', and no matter, giving the temperature increase. This is then 
simply related to the fractional increase in entropy density s'/s. 

4/3 4 [M\ ( M \ 2 ( T D \- x 

1 + - Y M I — ) = 1 + 4.87 ^— . 6 

3 \TJ V 10 PeV/ V 10 MeV/ yJ 

Here we gave the temperature in MeV to be suggestive. The temperature at decay 
must be larger than this so as not to disrupt the nucleosynthesis era. To be safe, 
To >10 MeV is warranted. 

To be specific, consider a gravitino with mass 1 keV, and thus VLz^h 2 ~ 1. Assum- 
ing we want Qh 2 = 0.15, we must dilute by a factor of s'/s = 6. For this model we 
find that the temperature increases by a factor of T'/T w 1.8 during the decay. We 
now find a simple relation between the mass of the messenger and the temperature 
at which it decays, 

2 (7) 
.10 PeV/ V 10 MeV/ V ; 

All that remains to be done in this model is to adjust the coupling between the 

messengers and the Standard Model particles to enforce this constraint on the lifetime 

r of the messengers, 



1 1 Af rf 45 / Tjj \ 9o T 

r =- = _ = — ►( —£— ) ^l.UxlO 29 



r 2H Tg\/27r 2 ^ VlOMeV/ 10 PeV 
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Note that here is the effective number of degrees of freedom at the decay, which 
for temperatures between 1 and 100 MeV is g* = 43/4. Therefore, the coupling of 
messenger particles to the standard model particles must be extremely suppressed. 
This is actually a desirable feature so as to not spoil the flavor-independence of the 
gauge-mediated supersymmetry breaking []nj. 

If the messengers decay directly to light particles, we might expect T ~ g 2 M, 
which would mean there is an extreme fine tuning of the coupling g. However, if the 
interaction is mediated by heavy particles of masses Mx, the lifetime is 

^ _g^M*_ ^ / T D A 2 ^ f_l_V f M \ 5 / M x \ ~ 4 

~ 1927T 3 M x \10 MeV/ ~ ' ^0.01 / UOPeW UO 12 GeW ' lJ 

Thus, a heavy particle of a mass of 10 12 GeV coupling the messenger sector to the 
standard model gives roughly the right lifetime. This is somewhat lower than the 
GUT scale, but still far above the messenger scale, so these interactions should not 
invalidate the relic density prediction. 



4 Entropy production: full evolution 

In the previous section we assumed that the messengers decay instantaneously. This 
approach allows us to study entropy production without lengthy calculations. The 



results from a more accurate calculation are similar, as we will show ]TJ|. In reality 
of course, the decay is exponential, and we should follow the evolution of matter and 
radiation through the decay process to accurately determine the increase in entropy, 
and the relative change in gravitino temperature. In this section we will do just that. 
We find a similar increase in entropy as for the simple calculation, though we also 
show that the temperature never increases, it merely decreases more slowly. 

We now write down the equations for the evolution of the radiation and matter 
densities, and the expansion scale factor. We assume that the universe is flat, with 
no cosmological constant, as is appropriate for early times. The equations are 

p M + 3Hp M = -Tp M , (10) 

/ 27r 2 \V3 

S 1/3 S = f—^J Yp M R\ (11) 

= ^(PM + Pn), (12) 
where S = sR 3 is the comoving entropy density. The first equation is solved trivially, 

PM®=PM&)(j^) 3 e- r(t ~ U \ (13) 
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and we solve the remaining two equations numerically. We recover the entropy in- 
crease from the solution, 

s> _ S(t f ) 



s S(U) ' 

The effective temperature decrease of the gravitino is also easily recovered as 

T R (t f )/T R (U) '--'^ 



T 3/ 2(tf)/T 3/2 (ti) 



(14) 



(15) 



in other words the radiation temperature decreases less than the gravitino tempera- 
ture. This indicates that for a given mass, the gravitino behaves more like cold dark 
matter than expected. 

Considering the results of the previous section, we choose two fiducial cases to run 
numerically. We use decay temperatures of 10 and 100 MeV, requiring M ~ 14, 45 
PeV respectively, and we set the decay times as appropriate. Evolving Eqs. ([10HT2|) 



for these parameters, we arrive at Fig. 1. We start the evolution at a temperature of 
10 GeV, and take g*(T) for entropy from Ref. |17]]. The onset of matter domination 
and the entropy-producing decay are both clear. By a temperature of a few MeV, the 
decay is all but complete, and nucleosynthesis can proceed as normal. We note that 
in principle, an epoch of matter domination such as is present in our model affects the 
evolution of the density fluctuations responsible for structure formation. In our case, 
the effects occur at wavelengths far too short to be observed in the cosmic microwave 
background or in large scale structure. 



5 Discussion 

We have constructed a model for gravitino warm dark matter that has the appropriate 
relic density. Entropy is produced by the decay of the messenger particles. The 
lifetime of the messengers is of the correct order if the decay happens through heavy 
particles of 10 12 GeV or so. 

For the keV gravitino we consider, the temperature, and thus the momenta, are 
decreased relative to the radiation by a factor of approximately two. This modifies the 
free-streaming scale Rf s of the gravitinos, and thus the effect on structure formation. 
This length scale is proportional to the velocity of the particles at the epoch of 
matter-radiation equality, and thus Rf s oc 23/2/777,3/2, normally given by JT3], [0| 

where in the second relation the expression for the relic density is substituted for 
mass. In our case the relic density is different, and we must be careful in replacing the 
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Figure 1: Evolution of entropy density. Here we show the decay of messengers at 
temperatures of 100 MeV (solid curve) and 10 MeV (dashed curve), producing an 
entropy increase of about a factor of eight. The onset of matter domination at a 
temperature of about 1 GeV (100 MeV) is clearly seen, and the decay at 100 MeV 
(10 MeV) is also clear. 
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expression. With entropy production, the gravitinos are colder at matter-radiation 
equality than a similar constituent would be without entropy production. This free- 
streaming scale comes about at the epoch of matter-radiation equality, where the 
temperature is a few eV. The gravitinos at this time are non-relativistic, so we can 
simply decrease their velocity by the decrease in temperature. We thus find that 
roughly speaking, 

This expression is identical to the previous one. If we use the quantities m 3 / 2 and 
^3/2 h 2 only, the expression for the free-streaming scale is unaffected by entropy pro- 
duction, at least in this approximation. 

A more accurate approach accounts for free-streaming prior to matter-radiation 



equality || |12| . We simply quote the results here. 

-1.15 



i4~0.65(ft m /i 2 )°' 15 Mpc. 



As the current bounds on the mass of a warm dark matter particle are in the 
0.75 — 1.0 keV range, we consider a gravitino slightly more massive than this, in the 
range 1.0 — 1.5 keV. The beneficial qualities of the warm dark matter scenario for 
the most part persist for particles in this mass range P|. The entropy production 
required is then 

S - = 7.6 ( J^VV^! V 1 ^11. (19) 
s UOOy V °- 15 / keV K ' 

We have shown that entropy production of this order is not a serious difficulty, thus 
the gravitino is again a viable candidate for warm dark matter. 
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